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We investigate tunneling properties of collective modes in the polar phase of a spin-1 spinor 
Bose-Einstein condensate. This spinor BEC state has two kinds of gapless modes, i.e., Bogoliubov 
mode and spin-wave. Within the framework of the mean-field theory at T = 0, we show that 
these Goldstone modes exhibit the perfect transmission in the low-energy limit. Their anomalous 
tunneling behaviors still hold in the presence of superflow, except in the critical current state. In the 
critical current state, while the tunneling of Bogoliubov mode is accompanied by finite reflection, the 
spin-wave still exhibit the perfect transmission, unless the strengths of a spin-dependent and spin- 
independent interactions take the same value. We discuss the relation between perfect transmission 
of spin-wave and underlying superfluidity through a comparison of wavefunctions of spin-wave and 
the condensate. 



I. INTRODUCTION 

The perfect transmission through a potential barrier, 
which is frequently referred to as the anomalous tun- 
neling phenomenon, is now widely recognized as a fun- 
damental property of the Bogoliubov excitation in a 
Bose superfluid [ll4l6||. This phenomenon is quite dif- 
ferent from the simple tunneling problem discussed in 
a one-particle quantum mechanics, where the transmis- 
sion probability always vanishes in the low-energy limit. 
The anomalous tunneling is known to also occur in the 
presence of superflow, except in the critical supercurrent 
state, where the tunneling of Bogoliubov mode is accom- 
panied by finite reflection. 

In Refs. 0, EMI, we have extended the previous 
work for the anomalous tunneling behavior of Bogoliubov 
mode to spin-wave excitations in the ferromagnetic phase 
of a spin-1 spinor Bose-Einstein condensate (BEC). We 
clarified that, not only the Bogoliubov mode, the gapless 
transverse spin-wave also exhibits the perfect transmis- 
sion in the low-energy limit. In the presence of a finite 
superflow, the perfect transmission of the spin-wave oc- 
curs when the magnitude of the spin-wave momentum 
k equals the supercurrent momentum q. On the other 
hand, the longitudinal spin-wave with a finite excitation 
gap does not show the anomalous tunneling behavior. 

In this paper, we further extend our previous stud- 
ies [l2l. [T^4T6j to the polar phase of a spin-1 spinor BEC. 
Recently, two of the authors [l5[ clarified that the spin- 
wave in this phase also exhibits the perfect transmission 
in the absence of superflow. In this paper, we examine 
how this anomalous tunneling behavior is affected by su- 
percurrent. Within the mean- field theory at T — 0, we 
show that, in contrast to the case of gapless ferromagnetic 
spin-wave, the polar spin-wave always tunnel through a 
barrier without reflection in the low- energy limit. Al- 
though this property is the same as the anomalous tun- 
neling behavior of Bogoliubov mode, we clarify that the 
perfect transmission of polar spin-wave occurs even in 
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FIG. 1: Model one-dimensional tunneling system. The sys- 
tem is in the polar state <&(x,t) — (0, &o(x, t), 0) with a fi- 
nite superflow in the positive z-direction. We examine the 
transmission and refection of a collective mode injected from 
x — —00 by the potential barrier V(x). 



the critical current state, except for the case in which 
the spin- dependent interaction strength equals the spin- 
independent interaction strength. This critical behavior 
is quite different from the Bogoliubov mode, which is ac- 
companied by finite reflection. 

This paper is organized as follows. In SecHH we explain 
the outline of our formulation to deal with the polar state 
of a spin-1 spinor BEC. In Sec lIIIl we examine the tun- 
neling properties of spin-wave excitations and effects of 
superflow in the polar phase. 



II. MEAN-FIELD THEORY FOR THE POLAR 
STATE OF SPIN-1 SPINOR BEC 

We consider the polar state of a spin-1 spinor BEC, 
and examine the tunneling of a collective mode through 
a potential barrier V{x) put around x — 0. Assuming 
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the uniformity of the barrier in the y- and z-direction, 
we treat this problem as a one-dimensional system, de- 
scribed by the Lagrangian density [TH-flU , 

C =ih&(x,t)dt${x,t) 

h 2 



^(x,t)(-—dl + V(x))Hx,t) 

C 2/ 



+ ^p'(x,t) + ^F 2 (x,t) 



Here. 



<f>o(x,t) 
*-i(x,t) 



(1) 



(2) 



is a spinor Bosc field with mass m, where t) 
(j = ±1,0) describe three magnetic sublevels. p(x,t) = 
$ t (x,t)$(x,t) and F = & (x,t)S$(x,t), respectively, 
describe the particle density and spin density, where 
S = (S x , S y , S z ) are the S = 1 spin matrices (where 
we take the spin quantization axis parallel to the z- 
axis). Co = 47r h 2 (ao + 2a2)/(3m) and ci = Ai:h 2 {a2 — 
ao)/(3m), respectively, describe a spin-independent and 
spin-dependent interactions [3], where as is an s-wave 
scattering length in the channels of total spin S = 0, 2. 

In Ref. [16j, we have explained the detailed mean- field 
theory for the model system described by Eq. ([1}, to 
examine the tunneling problem in the spin-1 ferromag- 
netic spinor BEC. Since the previous formulation is also 
applicable to the present polar state, we only present the 
outline of our formulation here. For more details, we refer 
to Ref. 0. 

The time-dependent Gross-Pitaevskii equation for the 
condensate wavefunction $(a;,i) is given by 
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d< f>(x,t) 



f h+(x,t) ^f- 
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Hx,t), (3) 



h 2 dl/(2m)+V(x)+c p(x,t), h±(x,t) 



where h{x, t) = 

h(x,t) ± c\F x , and F± — F x ± iF y . As usual, we 
study the tunneling problem in the stationary state, 
as schematically shown in FigfTJ In this case, setting 
$(x,t) = e~ v */ ?i <l(x), we obtain the GP equation for 
the spatial part $(x) as 
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^±( x ) = ft.(x) ± c\F z , where p(x) 
condensate density. 



V(x) + cop(x) and 
= &(x)&(x) is the 
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FIG. 2: (Color online) Calculated condensate wavefunction 
Qq{x) — \&o(x)\e lvo< ' x ' > around the potential barrier V(x) = 
2copo exp [— (a;/£) 2 ] in the critical current state, (a) Ampli- 
tude |"I > o(a;)j. (b) Phase ipa(x). In panel (b), we set <po(0) = 0. 
The inset shows the relation between the phase difference 95 
(which is given in panel (b) [2l[) and supercurrent momentum 
q, where £ = h/ \/rfiCoPo is the healing length. The red point 
at q = 0.0574^ _1 (= q c ) represents the critical current state. 



The polar state becomes the most stable state when the 
spin-dependent interaction is antiferromagnetic (c\ > 0) . 
In this case, the condensate wavefunction <t(x) has the 
form 11911 . 





$(a;) = ( * o (a0 




(5) 



which leads to the vanishing spin density F = 
$t(:r)S$(x) H- The GP equation for the non-vanishing 
$o(:c)-component is then given by 



h(x)* (x) = 0. 



(6) 



In the supercurrent state, since the superfluid properties 
far away from the barrier are the same as those in the 
uniform system, one finds p = cqPq + h 2 q 2 / (2m) (where q 
is the supercurrent momentum and po is the condensate 
density at x — ±00). We show in Fig. [5] the calculated 
spatial variation of the condensate wavefunction <&o(x) 
in the supercurrent state. 

A collective mode in the polar phase can be con- 
veniently calculated by considering fluctuations of the 
condensate wavefunction around the mean-field value in 
Eq. ([5]). In the time-dependent GP equation ([3]), set- 
ting $>(x,t) = exp(~ipt/h)(<t(x) + <j>(x,t)) (where <fi = 
(</>+i, 0o ; 4>-i ) describes fluctuations around the mean- 
field value $>o(x)), and retaining terms to 0(<p(x)), we 
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obtain 



ifi dfofo*) =h fa Cq)Mx) + c^)^), (7) 
.d<p±i(x,t) 



dt 



=h(x,ci)(f> ±1 (x) +c 1 <S> 2 {x)(j)* 1 (x), (8) 



where h(x, c ,i) = + c ,i|$o(^)| 2 - Since Eq.© for 
0o is decoupled from Eq.® for <j)±\, one may separately 
treat them. The former equation gives the Bogoliubov 
excitation. Indeed, setting (f>o(x,t) — u{ y x)e~ lEt l n — 
v*(x)e +lEt / h , we obtain the ordinary Bogoliubov equa- 
tion, 

fu{x)\ _ ( h{x, Co) ~c <I> 2 (x)\ fu(x)\ . . 
* ^(z) ) ~ [c^f(x) -h(x, c ) ) [v(x) ) ■ W 

On the other hand, one may interpret 4'±i{x,t) as fluc- 
tuations of S z — ±l-components around $±i = 0, 
so that they give spin-wave excitations. Substituting 
<f)±i{x,t) = <f) ±1 (x)e TlEt/h into Eq. ©, one obtains the 
"spin-wave equation," 
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4>+i{x) 



h(x,ci) 



ci^(x) 
-h(x,ci] 



0+1 0) 



(10) 



When (cj) + i(x),(j)-i(x)) is a solution of Eq. (ITUl) . 
(0_i(a;)e i£;t / R ,0 + i(x)e- l ' Bt / ?i ) also satisfies Eq. ©, 
which means that the spin-wave excitations in the po- 
lar state are doubly degenerate. 

The set of GP equation (J6]) and the Bogoliubov equa- 
tion §§§ is the same as that in the scalar BEC. Thus, 
the Bogoliubov excitation in the polar phase of a spin-1 
spinor BEC has the same tunneling properties as that 
in the scalar BEC. That is, the perfect transmission oc- 
curs in the long- wavelength limit, and this anomalous 
tunneling phenomenon still holds in the presence of a fi- 
nite superflow, except in the critical current state. In 
the critical current state, the transmission of low-energy 
Bogoliubov mode is accompanied by finite reflection. 

The spin- wave equation (I10[) reduces to the Bogoliubov 
equation ([9]) when c\ = cq [22j . When this relation is sat- 
isfied, spin- wave excitations are found to also exhibit the 
anomalous tunneling behavior as in the case of Bogoli- 
ubov mode. However, one cannot use this analogy when 
ci 7^ Co. The goal of the next section is to clarify what 
happens in this case. 

In considering the supercurrent state, we need to 
choose the supercurrent momentum q so that the Landau 
instability will not occur. In the uniform system, Eqs.([9]) 
and (|10p , respectively, give the Bogoliubov and spin- wave 
excitation spectra as, 
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Thus, the stability condition for the supercurrent state, 
which is guaranteed when Eqs. (fTTj) and (IT^t are always 
positive, is obtained as 



< Min[l, y/ejoo] 



(13) 



For later convenience, we introduce the dimensionless 
variables, E = E/(cqPq), x = x/^, k — q = q£, 
V = V/{c p ), = ^i/y/p^, ~$ t = <kU/ Pa, i = tcop /h 
and c = ci/co, where £ = h/y/mcopo is the healing 
length. Using these, we examine the low-energy tunnel- 
ing properties of the polar spin- wave in the next section. 
For simplicity, we omit the bar in what follows. 



III. ANOMALOUS TUNNELING OF 
SPIN- WAVE MODES IN THE CURRENT 
CARRYING POLAR PHASE 



We numerically solve the spin- wave equation (jlOp . to- 
gether with the GP equation ([6]), in the presence of a 
finite superflow with the momentum q. As shown in 
Fig [1] we consider the case that an incident spin- wave 
with the energy E is injected from x <C — 1, and tunneling 
through the barrier V(x) around x = 0. The wavefunc- 
tion (</>+i (x), 4>^i{x)) obeys the boundary conditions, 



+i i x \ -j. (ah e +lqx \ i 



A 



-\-iqx 



a k3 e 
Pk 3 e- iqx 



e lfe ^(x« -1), 



(14) 



hi] _ T (a kl e +lqx \ ikix 
f-J \^e- iqx ) 

+ Bh^Z)e ikiX (x^+l). (15) 

Here, each term in Eqs. (fT4"|) and (fT5")) is an independent 
solution of Eq. (fTUf in the absence of barrier V{x). That 
is, the coefficients a k and (3 k are given by 



Ok 

Pk 



1 



E-c-qk~k 2 /2\ 2 
c 



E 



qk - k 2 /2 



(16) 
(17) 
il 2 - 

1 following the Bogoliubov mode [23j. The 
momenta kj (J = 1 ~ 4) are the solutions of 

k 4 + 4(c - q 2 )k 2 + 8kqE - 4E 2 = 0. (18) 

Equation (fT8|) has two propagating solutions (j = 1,2), 
as well as two damping solutions (j = 3,4). In the low 
energy region, they are given by 



where the normalization condition is imposed as \(t>+i 

l<Ml 2 



E 



k 2 



Vc + q 
-E 



(19) 
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Once the coefficients (T, R, A, B) are determined, the 
transmission probability can be calculated in the same 
manner as in the ordinary tunneling problem. From 
Eq. ((5J, we find that the generalized density n a (x,t) = 
\<fi±l(x, t)\ 2 — \(f>—i(x, t)\ 2 satisfies the continuity equation 
dtn s + d x J s = 0. Thus, in the stationary state, the flux 
density 

J s = -lm[cp* +1 d x 4> +1 + 0-10^*1] (21) 
m 

is a conserving quantity. Using this, one may conve- 
niently introduce the transmission probability r, as well 
as the reflection probability r, as 



If, 
w(k 2 ) 



\R\ 



(22) 
(23) 



where w(k) = k[\a k \ 2 + \f3 k \ 2 ] + q[\a k \ 2 - \fa\ 2 ]. J> and 
Js describe the flux densities of incident and reflected 
waves, respectively, that are calculated from the first and 
second terms in the RHS of Eq. (fT4)) . respectively. The 
flux density J* of the transmission wave is obtained from 
the first term in the RHS of Eq. (fl5|l. 

Figure |3] shows the transmission probability r of the 
spin-wave through the potential barrier, 



V(x) = Voe-<-* 



A) 



(24) 



In this figure, we find that the anomalous tunneling (per- 
fect transmission) occurs in the supercurrent state of the 
polar phase, irrespective of the value of c = Cx/cq. The 
perfect transmission can be seen even in the critical cur- 
rent state, except when c — 1, which is quite different 
from the case of Bogoliubov mode, where the tunneling 
is accompanied by finite reflection (r > 0). 

We briefly note that the tunneling properties of the 
spin- wave mode shown in Fig|3] are quite different from 
the case of ferromagnetic spinor BEC . In the latter 
case, the gapless transverse spin-wave mode exhibits the 
perfect transmission in the low-momentum limit only in 
the absence of supercurrent. In the supercurrent state, 
the perfect transmission occurs at the momentum whose 
magnitude equals that of the supercurrent momentum. 
The other spin-wave mode with a finite excitation gap 
does not show the anomalous tunneling behavior. 

As mentioned previously, the spin-wave equation (|10[) 
reduces to the ordinary Bogoliubov equation when c = 
c i/co = 1- Thus, the absence of the perfect transmission 
in the critical current state shown in Figj3] (b) is due to 
the same mechanism as that in the case of Bogoliubov 
mode. In the latter case, in addition to the zero-energy 
phase mode solution given by 
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v{x) 



*S0»0 



(25) 
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FIG. 3: Transmission probability r of the spin- wave in the 
polar phase, as a function of incident momentum k. We take 
Vo — 2copo and A/£ = 1, that give the same potential barrier 
as that used in Fig. [2] 



the Bogoliubov equation (fSJ has the other zero-energy 
solution in the critical current state [Toj . 



u(x) 
v(x) 



(26) 



where ip is the phase difference of the condensate wave- 
function <&o(x) between x = ±oo [2l[. This solution 
physically describes density fluctuations near the barrier, 
leading to the suppression of perfect transmission [To| . 
Similarly, Eq. (|10[) always has the zero-energy solution, 



<t> +x (x) 
6-i(x) 



*o(a0 



(27) 
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which describes the Goldstone mode associated with the 
broken spin rotational symmetry in the polar phase. In 
addition to this, one also has the other zero-energy solu- 
tion in the critical current state, given by 

( 4>+i{x) \_( d v $ (x) \ . . 

While spin fluctuations described by SF± = \/2~(3>o<A±i + 
<J>50 T i) vanishes in the case of Goldstone mode in 
Ed. ([27]), Ea. (155)1 gives SF± ^ 0. Thus, in contrast to 
the case of Bogoliubov mode (where density fluctuations 
destroy the perfect transmission in the critical current 
state), the absence of the anomalous tunneling of the 
spin-wave when c = 1 is due to the emergence of spin 
fluctuations. 

We note that Eq. ((55)) does not satisfy Eq. (fTU)) unless 
c = 1. Because of this, when c ^ 1, the spin- wave can 
tunnel through the barrier without being disturbed by 
spin fluctuations, even in the critical current state, as 
shown in Figs. [3] (a) and (c). 

We also note that, as pointed out in Ref. ||24|], the emer- 
gence of the zero-energy density mode in Eq. (|5o]) is 
deeply related to the instability of the superfluid state 
in the case of scalar BEC. In the polar state of a spin-I 
spinor BEC, while both low-energy density fluctuations 
and spin fluctuations exist when c = 1, the former fluc- 
tuations only appear when c ^ 1 in the critical current 
state. This implies that the mechanism of supercurrent 
instability above the critical current state might be some- 
how different between the cases of c = 1 and c ^ 1. 

Equation (|57[) indicates that the general property of 
the zero-energy Bogoliubov mode, that each component 
of the wavefunction (u(x),v(x)) coincides with the con- 
densate wavefunction 3>o(x) 0] (See Eq.([25)l). also holds 
in the spin-wave. In this regard, it has been also pointed 
out in the scalar BEC that such a coincidence also ex- 
ists, when the Bogoliubov excitation has a small but fi- 
nite momentum (E cx p > 0) [7] . To see if this property, 
which we call "supercurrent behavior" , is also applicable 
to the spin-wave, we consider the tunneling problem for 
a (5-functional potential barrier V(x) = VqS(x). In this 
simple model, solving the GP equation © analytically, 
one obtains the condensate wavefunction in the super- 
current state with momentum q as [H, @] , 

$ (x) =e^ x+e «\ 1 <(x)+iq] 

=* < (x,< Z )(x<0), (29) 

$ (x) =e i ^- 6 ^[ 1 >{x)-iq] 

=*>(*,?)(* >0), (30) 

where J^~(x) is J q (x) for x < 0, J^(x) is ~/ q (x) for x > 0, 
and exp(i0 q ) = [j q (0) - iq}/ JjqHP) + 1 2 - Herc > lq( x ) = 

y/l - q 2 tanh (y/l - q 2 (\ x\ +xq)), and Xq is determined 
so as to satisfy d x $ (+0) - d x $o(-Q) = 2Vb$ (0) 0. In 
the low momentum region (|p|, \q\ <C 1), one obtains the 



following form: 

(<t>+l \„ (K^<{X-,p + q)\ , . ( ) 

{r-J-{^ < (X + ,p~q))^ X<Uh 

with the accuracy of 0(jp|), and 0(|g|). Here A± are 
unimportant numerical factors which depend on p and q; 
X± are functions of x. We note that the relations (|3"T)) 
and (|32[) cannot be extended to q — q c when c = I. If 
(I3T1) and ((32]) held for q = q c when c = 1, it would lead 
to partial transmission of spin-wave at low-energy limit. 

"Supercurrent behavior" was first proposed as an ex- 
planation of perfect transmission of the Bogoliubov mode 
through the potential barrier Q . This property was suc- 
cessfully proved in the transverse spin- wave of the ferro- 
magnetic spinor BEC for an arbitrary shape of the bar- 
rier [l6| . The Bogoliubov mode and the spin- wave in the 
present polar phase certainly show "supercurrent behav- 
ior" in the presence of the (5-function potential barrier, 
but it is fair to say that the applicability of "supercur- 
rent behavior" as an explanation of perfect transmission 
of these two modes through an arbitrary shape of the bar- 
rier is currently restrictive compared with the transverse 
spin- wave in the ferromagnetic spin-I BEC. 

IV. SUMMARY 

To summarize, we have investigated tunneling proper- 
ties of collective excitations in the polar phase of spin- 
I spinor BEC. We showed that the low-energy spin- 
wave can tunnel through a barrier without reflection, 
even in the presence of a finite superflow. In the crit- 
ical current state, tunneling properties depend on the 
ratio c = Cx/co, where Co and c\ are the strengths of 
a spin-independent and spin-dependent interactions, re- 
spectively. When c — 1, the perfect transmission does 
not occur in the critical current state. Otherwise, the 
anomalous tunneling behavior continues to exist in the 
critical current state. 

In the case of Bogoliubov mode, the anomalous tun- 
neling phenomenon does not occur in the critical current 
state. In the critical current state, a density fluctuation 
mode appears around the barrier, which suppresses the 
perfect transmission of Bogoliubov mode. In the spin- 
wave case with c = I, spin fluctuations appear around 
the barrier in the critical current state, leading to the 
breakdown of the perfect transmission of the spin-wave. 
When c ^ 1, such spin fluctuations do not appear, so 
that one obtains the anomalous tunneling phenomenon 
of the spin-wave in the critical current state. 

We note that the anomalous tunneling behavior of the 
spin-wave mode in the polar phase is somehow different 
from the case of ferromagnetic spinor state. In the latter 
case, the gapless transverse spin- wave excitation exhibits 
the perfect transmission, not in the low-momentum limit, 
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but in the case when the magnitude of the spin- wave mo- 
mentum k equals that of supercurrent momentum q. The 
longitudinal spin-wave with a finite excitation gap does 
not show the anomalous tunneling phenomenon. In this 
case, the perfect reflection occurs in the low-momentum 
limit, as in the case of the ordinary tunneling of a quan- 
tum mechanical particle. 

However, when the perfect transmission occurs, the 
wavefunctions of Bogoliubov mode, polar spin-wave, and 
ferromagnetic transverse spin-wave, all have the same 
form as the condensate wavefunction. Through a se- 
ries of papers and the present work, we have shown that 
anomalous tunneling is a common property of all gapless 
collective modes in scalar and the spin-1 spinor BECs. 
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Appendix A: Derivation of Eqs. (pTTjl and (f32j) 

Expanding the condensate wavefunction in Eqs. (|29[) 
and (|30p to 0(q), we obtain 



where {X,Y} = Xd x Y - [d x X)Y, and 



*>(x,?) ^e 1 



J K (x) - iq 



J > (x) + iq 



7 <(x)-j(0) 



7 (0) 
7 >(z)- 7 (0) 



7 (0) 



(x < 0), 

(Al) 

0= > o), 

(A2) 



where r y < (x) is j(x) for x < 0, 7 > (a;) is 7(2:) for x > 0. 
Here, 7 (x) = tanh(|a;| + xq). In this limiting case, the 
boundary condition at x = (d x &o(+0) — d x $o(—0) = 
2Vo$ (0)) gives V = [1 - 7 2 (0)]/ 7 (0). The goal of this 
appendix is to show that the components 4>±\{x) in the 
wavefunction of the spin- wave with \p\ -C 1 reduces to 
the forms (gl]) and (fX2|) . 

We introduce the functions S(x) and G(x), given by 

S(x) =<j) +1 {x)e- lipa{x) - (f>*_ 1 {x)e +lVa{x \ (A3) 
G(x) =<j) +1 {x)e- lipa{x) + (f>t 1 {x)e +lVa{x \ (A4) 

where <fio(x) is the phase of the condensate wavefunction 
$ ( x ) = |$ (a;)|e^ ( x ). In the case of 5-functional bar- 
rier [3j, we obtain tpo{x) = qx + sgn(a;)[tan~ 1 (7 9 (x)/g) — 
tan _1 (7 g (0)/g)]. Equations for S(x) and G{x) are ob- 
tained from Eq. ([TUl) as 



EG{x) =H(x)S{x) - iq 



{\* Q {x)\,G(x)} 



|$ (x)| 3 : 
ES{x) = [H(x) + 2c|$ (a;)| 2 ] G{x) - iq 



(A5) 
{\$ (x)\,S(x)} 



H 



w + 2>w +y ^- 1+ y + i $0 ^i 3 - (A7) 



The boundary condition in terms of S(x) and G(x) is 
given by 



fc i 1 cikix j_ R I ^ k2 \ 
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ki 



Gk 2 



+ A (g*^e ik3X (x «-l), (A8) 



B I J** ) e 4fc ^(x hi !. 



(A9) 



where (§ k ,G k ) = (1 + k 2 /(4:c),(E - qk)/(2c)) for fc = 
fe ll2l and (5fc,G fc ) = ((£? - gfc)/(2c), fc 2 /(4c)) for k = 
^3,4 [111- I n the low energy regime (£" <C 1) with small 
supercurrent momentum \q\ <C 1, one may approximate 
the spin- wave momentum to k\ ~ E / ^/c (= p). Then, in 
expanding coefficients in terms of p and q as 



U'=0 



/ T (M')\ 



(A10) 



the boundary condition of the lowest order of p and q is 
given by 



1 + \ 

A (0,0) e V5c ) < ( AU ) 

71(0,0) \ 

B (0,0) e -2^J (*»+!)■ ( Al2 ) 



On the other hand, expanding S 1 and G in terms of p 
and g as 



00 



(A13) 
(A14) 



Z,Z'=0 



(A6) 

we obtain equations for the coefficients in Eqs. (|A13l) and 
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(1AT4|) as, 

h {x)S^- 0) {x) =0, 
[h (x)+2 n 2 (x)} G^ Q \x) =0, 

[h (x) + 2c 7 2 (x)] G^°\x) =^S^°\x), 
ho(x)S<M(x) 



[h (x)+2 cl 2 (x)] G (0 < 1} (x) =j 



. { 7 (x),G( - >(x)} 
7 3 (x) 

.{ 7 (x),S( < )(x)} 
7 3 (x) 



(A15) 
(A16) 
(A17) 
(A18) 

(A19) 
(A20) 



where h (x) = -(l/2)o£ + V S(aO - 1 + 7 2 (x). 

Equation ()Al5j) has two independent solutions f c (x) = 
j(x) and f (x) = xj(x) + sgn(x)[7(x) — 7(0)]/7(0), both 
of which satisfy the boundary condition at x = 0. As 
a result, S(°'°)(s) = C^f e (x) + D<- ^f o (x) follows, 
where G^ 0,0 ) and Z)( ' ) are coefficients. Taking |x| — >■ 
oo and considering the boundary conditions (|A1 1|) and 
(TAT21) . we obtain Geo* ) = T(°>°) = and = 

0. 

To solve Eq. (|Af 6[) , it is convenient to rewrite this equa- 
tion in the form 



(1-7 2 ) 
2(1 + 2c) 



dj 2 
4c 



l- 7 2 



2 dG(°'°) 
c?7 

G(°'°) = 0, 



(A21) 



where we have used the relation d x j(x) = sgn(x)[l — 
7 2 (x)]. The solution of (|AT6)) is given by G (0 ' \x) = 
aPp(-f(x)) + bQ 1 l(j(x)), where Pf} and Q2 are the associ- 
ated Legendre functions [26[, with v = (—1+^/9 + 16c)/2 
and 77 = 2y / c. When 77 is a non-integer, the asymptotic 
behaviors of the associated Legendre functions are given 
by 



KM*)) 



T(l-r,) 



2 tan(r?7r) T(l - rj) 



(A22) 
(A23) 



When we take b = — a(2/ir) tan(?77r)(= ab), the func- 
tion GQ'°\x) = a[P2(rr(x)) + 6Q3(7(a))] 

converges. 

However, imposing the boundary condition at x = 0, 
d x G^(+0) - 9 x G(°>°)(-0) = 2VbG(°<°)(0), one obtain 
a = 0, leading to G(° ,0 )(x) = 0. When r\ is an integer, 
while PJ( , y(x — > ±00)) converges as exp (— r)Xa)T(v + 
?7 + l)/[r)T{v -77 + I)] exp (-r)\x\), Q2("f(x)) diverges as 
oc exp (r]\x\) (\x\ S> 1). Thus, we need to take b = 0. 
However, the boundary condition at x = again gives 
G(°<°)(x) = Thus one concludes G<-°^(x) = 0. Con- 
sidering (TATTt and (|Al~2|) . we have A' ' ) = = 0. 
The boundary conditions (| A8|) and (| A9|) are then given 



by 

5 \ _ / T (o,o) + p[jR (i.o) + ia .^ _ ^(o.o))] + gjR (o,i)> 
Gj - ^ p/(2^)+^ 1 )e 2 ^ 

(<b«-1) 
(A24) 

S \ _ Mo,o) +p(T (i,o) +ttT (o,o) ) + (;r (oi)N 

Gj - ^ P /(2^E)+qB^' 1 h~ 2 ^ cx 

(x > +1) 
(A25) 

within the accuracy of 0(p) and 0(q), where we took 
p\x\ < 1. 

Noting that G (0 ' 0) (x) = 0, one finds that Eqs. (|AT7|) 
and (|A19|) are the same as Eq. (IA15|) . Thus, S^ 1,0 )^) 
and S( '^(x) are both given by the linear combination of 
fc(x) and fo(x). Using this, we may write S(x) in the 
form, within the accuracy of 0(p) and 0(q), 



S(x) =TC°.°)/ B (x) + p (g^ >/ o (x) + D^f a (x)) 

+ q[c iQ ^f c (x)+D^f (x)). (A26) 



G< M ') and DM") (1,1' = 0,1) are coefficients. 

Comparing Eqs. (|A26[) for |x| ^> +1 with Eqs. 
(IA241) andfTMSjl^ehave^C^^t ' ),^^' ),^^' 1 )) = 
(1, 0, z, 0), G(^°) = TC 1 ' ) = and T^°) - T^ 1 ) = 

- i?(°^)) = + 1/7(0)]. Using this result, 
one finds, 

S(x) =[1 + i{pC {lfi) + qCW) + ipx] 1 (x) 

+ sgn(x)ip[ 7 (x) - 7(0)]/7(0). (A27) 

As a result, we obtain, within the accuracy of 0(p) and 
O(q), 

S(x) =e ia ™y < (x, q)(x < 0), (A28) 
S(x) =e ia "<«*>(x, q)(x > 0), (A29) 

where ot M = pC^ + qC^ ^. 

Next, we evaluate G^ '^ and G (1 ^ 0) . Substituting 
5(0,0) _ intQ Eqs _ (TJ2Q)) and (|AT8)) . we obtain 

[/i (x) + 2c 7 2 (x)]G (0 ' 1) (x) =0, (A30) 

[/i (z) + 2c 7 2 (x)]G( 1 <°)(x) =Vc~ 7 (x). (A31) 

As in the case of G (0 ' 0) (which obeys Eq. (fAl6| ), Eq . 
(lA30)l gives G (0,1) = 0. Considering (TAT! and (TA251) . 
we have A^ - 1 ) = S^ ' 1 ) = 0. For Eq.JMEJ, we have, 
within the accuracy of 0(p) and 0(q), 

G(x) =p[a^P?(j(x)) +b^Ql{ 1 {x)) 

+ G^)(x)], (A32) 

where 

sgn(x) 



Gi l ' 0) («) 



2^fc 



ll{l{x))l{x)dx 



■\< l X S -^P-P?(7(x)) 7 (x)dx 



Q n Ai{x)). 

(A33) 
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X\ and X2 are constants and A is defined as [2] 



A = 2 2Tl e 2T > 7 



i r((i/ + g + l)/2)r((i/ + q)/2 + l) 
T((p-ri + l)/2)T((p-ri)/2 + iy 



(A34) 



Considering the asymptotic behavior of PJ('y(x)) and 

Q2(l( x )) f° r x ~ we find that G^ 1,0 ' (a;) converges. 

Since G(x) should converge for x = ±00 in this tun- 
neling problem, we have &T ) = a^ 1,0 ^b for 77 being a 
non-integer and &(L0) = for rj being an integer as in 
the same discussion for G^ 0,0 ^(x). We still have three 
unknown parameters (a^ 1,0 -*, x\, x%), but these three pa- 
rameters are determined from three boundary conditions: 
First one is the boundary condition at x — given by 
d x G{+0) - d x G(-0) = 2V G(0). Second and third ones 
are those for x = ±00, given by G^'°\x — ±00) — 



l/(2y / c) where we used (|A24[) and (|A25|) . 

From Eqs. (MJ), flM}, (El, (IM) . (TM9]),and 

(|A32j) . we have 



<z)/2\ 
?)/2 J 



(x<0), 
(A35) 

/+e^.^ > (X +) p + g )/2A 

^ -^^(X-.p - q)/2 J ' lX " UJ ' 

(A36) 



within the accuracy of O(p) and O(q) if we use X± = 
x±G(x)/[l - 7 2 (x)]. Taking A± = ±e ia ^"/2, we obtain 
Eqs. ([3T ]) and ([32 ]) . 
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